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Abstract 

We present a semiclassical three-fluid model for a Bose-condensed mixture 
of interacting Bose and Fermi gases confined in harmonic traps at finite tem- 
perature. The model is used to characterize the experimentally relevant be- 
haviour of the equilibrium density profile of the fermions with varying compo- 
sition and temperature across the onset of degeneracy, for coupling strengths 
relevant to a mixture of 39 K and 40 K atoms. 
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The achievement of Bose-Einstein condensation (BEC) in trapped gases of 87 Rb 
[1], 23 Na [2] and 7 Li [3] has provided new impulse to the study of many-body and 
quantum statistical effects in dilute fluids at very low temperature. The formation 
of coexisting condensates by sympathetic cooling in a mixture of Rb atoms in two 
different internal states has also been achieved [4] . Trapping of fermionic species has 
been reported for 6 Li [5] and 40 K [6]. Trapped mixtures of bosonic and fermionic 
species are expected to become accessible to experiment in the near future. 

The density profiles of the separate fermionic and bosonic species in such mix- 
tures can in principle be experimentally resolved. With this perspective we present 
in this work a semiclassical three-fluid model extending our earlier studies of the 
condensate fraction and internal energy of a trapped interacting Bose gas [7, 8] to 
a Bose-condensed mixture of interacting Bose and Fermi gases confined in a spher- 
ically symmetric trap at finite temperature. We assume that the mixture is in full 
thermal equilibrium: the efficiency of a Bose condensate in cooling slow impurity 
atoms has been discussed by Timmcrmans and Cote [9] . The main emphasis of our 
calculations is on the evaluation of the density profile of the fermionic component 
with varying temperature and composition for values of the coupling strengths rel- 
evant to a mixture of 39 K and 40 K atoms. We also identify two parameters which 
govern the behaviour of the profile. In the limit of the Thomas-Fermi approxima- 
tion at zero temperature we recover the equations used by M0lmer [10] to describe 
the ground state of mixtures of bosons and fermions. Studies of the thermodynamic 
properties of ideal Fermi gases in harmonic traps have been reported by Butts and 
Rokhsar [11] and by Schneider and Wallis [12]. 

Interaction effects are very small in the normal phase but become significant 
as condensation induces a density increase of the bosonic component near the cen- 
tre of the trap. For values of the (repulsive) coupling strengths in the range of 
present interest, the condensate pushes both the bosonic noncondensate and the 
fermionic fluid towards the periphery of the trap. With regard to the latter fluid, 
such 'squeezing' drives an increase in its chemical potential. As we shall sec below, 
these effects on the fermionic density profile are especially striking in boson-rich 
mixtures, but rather low temperatures relative to the BEC transition need to be 
reached. Of course, an increase of the fermion-boson coupling, as may be achieved 
by modifying the scattering lengths with external fields [13, 14] and by changing 
the parameters of the trap or the components of the mixture, will enhance the 
BEC-induccd changes in the fermionic fluid. 

In the following we shall assume that the number of bosons in the trap is large 
enough that the kinetic energy term in the Gross-Pitaevskii equation for the wave 
function ^(r) of the condensate can be neglected (see e.g. [15]). This corresponds 
to the so-called Thomas-Fermi limit and is in general a good approximation except 
in the immediate neighbourhood of the BEC transition temperature. It yields the 
strong-coupling result 

* 2 (r) = [w, - V b ext (r) - 2gn nc {r) - fn f (r)]/g (1) 

when the quantity in the square bracket is positive, and 5 ,2 (r) = otherwise. Here, 
g = Aith ab/rrib and / = 2ttH 2 a / / 'm r with a/ and the boson-boson and boson- 
fermion s-wave scattering lengths and m r — m^mf / {mi, + m,f) with mi, and rrif 
the atomic masses; /if>(T) is the chemical potential of the bosons at temperature 
T, Vfr Xt (r) — nibLd b 2 r 2 /2 is a spherically symmetric external potential confining the 
bosons, n nc (r) is the average distribution of non-condensed bosons and rif(r) is 
that of the fermions. The factor 2 in the third term on the RHS of eqn (1) arises 
from exchange [16] and we have neglected a term involving the off-diagonal density 
of non-condensed bosons. 
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As already proposed in early work on the confined Bose fluid [17 - 19], we treat 
both the non-condensed bosons and the fermions as ideal gases in effective potentials 
V b ^ (r) and VJ^ (r) involving the relevant interactions. We write 



V b eff (r) = V b ext (r) + 2g* 2 (r) + 2gn nc (r) + fn f (r) 



and 



Vf s {r) = Vf xt {r) + /* 2 (r) + fn nc (r) 



(2) 



(3) 



with Vf Xt (r) — rrifUj'jr 2 /2. We are taking the fermionic component as a dilute, spin- 
polarized Fermi gas: the fermion-fermion interactions are then associated at leading 
order with p-wave scattering and are demonstrably negligible at the temperatures of 
present interest [11, 20]. We may then evaluate the thermal averages with standard 
Bose-Einstein and Fermi-Dirac distributions, taking the non-condensed particles to 
be in thermal equilibrium with the condensate at the same chemical potential ^ b {T) 
and the fermions at chemical potential fj,f(T). In the semiclassical approximation 
we obtain 



n f( r ) = ^3 / d3 P { ex P 



p- 



2m/ 



+ V? ff (r)- M ) /k B T 



1 



(4) 



and 



n nc {r) = ^ J d 3 p jexp 



2m b 



+ V b eff (r)- f i b ) /k B T 



- 1 



f2nm b k B T\ 3/2 ^ expj-j (V b eff (r) - n b )/k B T] 

' .7>1 



h 2 



J 



;3/2 



(5) 



The chemical potentials are determined from the total numbers of bosons and 
fermions, 



I 



d 3 r [* 2 (r) +n nc (r)} 



and 



N f 



I* 



■n/(r) 



(6) 



(7) 



These equations complete the self-consistent closure of the model. 

Before presenting relevant illustrative examples of the full numerical solution of 
the set of equations (1) - (7), it is useful to discuss a simplified form of the present 
three-fluid model (see also [7]). This is obtained by introducing an approximation 
which preserves only the repulsions exerted by the condensate: namely, we set to 
zero the last two terms in the RHS of eqns (1) and (2) and the last term in the 
RHS of eqn (3). Evidently, this approximation rests on the fact that both the 
non-condensed and the fermion component are very dilute gases and becomes all 
the more accurate as one moves towards boson-rich mixtures well below the BEC 
transition temperature. In this approximation we can (i) introduce a temperature- 
dependent scale length R — (2^i b /m b uJ b ) 1 ^ 2 defining the radius outside which the 
condensate density vanishes, and (ii) scale the effective potential acting on the 
fermions by writing it in the form Vj^ (r) = TiWfV^ \x) where x ~r/R and 



(x) 



i 7 [A+(1-A)a; 2 



b x2 



for 
for 



x < 1 

X > 1 



with 



gmjLjj 



(8) 



(9) 
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and 



2/j,),mfU)f 




1/2] 2 / 5 



(10) 



7 = 




= l5N c a b 



m b u> b 



Evidently, the dimensionless constant A controls the shape of the effective poten- 
tial seen by the fermions and hence their density profile: it depends only on (i) the 
ratio of the two relevant scattering lengths and (ii) the ratio of the spring constants 
of the two traps. Instead the parameter 7, which depends on the number N C (T) 
of bosons in the condensate through their chemical potential, controls the depth 
of the effective potential (in units oihujf). It is easily seen from eqn (8) that for 
A > 1 Vfff(x) decreases from the value 7A/2 at the centre of the trap towards its 
minimum value 7/2 at r — R, and increases thereafter because of the confinement. 
In this regime the fermions are squeezed away from the centre of the trap into a 
shell overlying both the outer part of the condensate and the non-condensate. On 
the other hand, for A < 1 the minimum value of Vj^(x) is 7A/2 at the centre of 
the trap: namely, in this regime the condensate merely raises the bottom of the 
confining well for the fermions without changing the sign of its central curvature. 

The calculations that we report below refer to the fermion density profiles in a 
trap with u)f = uJb = 100 s _1 . We shall first consider (left panels in Figures 1 and 2) 
the case a b — 80 and a/ ~ 46 Bohr radii for the boson-boson and boson- fermion ,s- 
wave scattering lengths, corresponding to the 39 K- 40 K mixture [21]. With the above 
values we find / ~ 0.57g and hence A < 1. We have also examined the regime A > 1 
by assuming / = 2g (right panels in the Figures). We have studied the dependence 
of the profiles on the composition of the mixture by considering two cases, namely 
(i) N f = 10 3 and N b = 10 6 (top panels in the Figures) and (ii) N f = N b = 10 4 
(bottom panels). The two characteristic temperatures of the mixture shown in 
the Figures are the BEC transition temperature ksT c = huj b (N b /^(3)) 1 ^ 3 for an 
harmonically confined ideal Bose gas and the Fermi temperature Tf. The latter has 
been calculated from the chemical potential of the fermions in the simplified model 
at zero temperature. The increase of Tf due to the boson-fermion interactions, 
relative to the confined ideal-gas value fc^T/ = hujf(6Nf) 1 ^ 3 , is quite large at boson- 
rich compositions (about 40% in Figure l.a and 86% in Figure l.c). 

Figure 1 compares our results for the density profile of the fermions at two 
temperatures below T c with those for a confined Fermi gas in the absence of all 
interactions. In this range of parameters the results obtained in the simplified model 
leading to eqn (8) are practically indistinguishable from those obtained from the full 
numerical solution of eqns (1) - (7). The interactions induce major distortions of 
the fermion density profile near the centre of the trap in the 39 K- 40 K mixture at low 
temperatures and at boson-rich compositions (see the curve in Figure l.a referring 
to T/T c = 0.1 and Nf = 10~ 3 7V 6 ). The transition from the first to the second 
regime of coupling strength A is clearly shown by the comparison between the left 
and right-hand panels in Figure 1. For boson-rich compositions at low temperature, 
the fermions are in fact almost wholly expelled from the central region of the trap 
(see Figure l.c). This agrees with the results of the ground-state calculations of 
M0lmer [10]. The density profiles of the bosons, that we do not show, are instead 
hardly affected by the presence of the fermions in the range of parameters that we 
have explored. 

It is also interesting to point out that for T = 0.1T C the conditions of Figures 
l.a and l.c correspond to temperatures of the same order as the Fermi temperature 
Tf. In this case the effects of the interactions are more important than those due 
to the quantum degeneracy of the Fermi gas. On the contrary, the conditions of 
Figures l.b and l.d at T = 0.1T C correspond to temperatures much smaller than Tf 
and the effects of quantum degeneracy are much more important, as is illustrated 
by the comparison with the density profile for an ideal Boltzmann gas at the same 
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temperature (dots in Figures l.a and l.b). 

Figure 2 shows the behaviour of the density nf(r = 0) of fermions at the centre 
of the trap as a function of temperature, for the same coupling strengths and com- 
positions as in Figure 1. In Figure 2. a we have reported in an inset the behaviour 
of the temperature derivative of n/(r = 0), to display its upturn occurring as quan- 
tum degeneracy develops at sufficiently low values of T/Tf. Major effects should 
be expected in boson-rich mixtures if the coupling strength can be driven into the 
A > 1 regime, as is seen from Figure 2.c. 

In summary, we have studied the equilibrium density profile of a fermionic fluid 
in a confined, Bose-condensed mixture of bosons and fermions in dependence of 
composition and temperature. The fermion density at the centre of the trap can be 
used to detect the onset of degeneracy in the Fermi gas with decreasing temperature. 
Although we have assumed spherically symmetric traps, the model can easily be 
extended to asymmetric confinements through a simple change of variables. 

As a final remark we notice that the density profile that we have calculated for 
the fermionic component immediately yields a semiclassical density of states p f (E) 
for the calculation of thermodynamic properties through the relation 

Pf (E) = 27r{2m f /h 2 ) 3 / 2 J d 3 r[E - Vf^r)] 1 ' 2 . (11) 

A similar relation holds, of course, for the density of states of the non-condensate 
[17-19, 7]. In the simplified model leading to eqn (8), that we have seen to be 
quite accurate compared with a fully self-consistent treatment of the mixture in 
the range of parameters that we have explored, the integral in eqn (11) is easily 
evaluated analytically. Calculations of thermodynamic properties for boson-fermion 
mixtures will be reported elsewhere. 
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Figure captions 

• Figure 1. Fermion density profile n/ (r) (in units of Og 3 , with a = (h/rribuj) 1 / 2 
and u the characteristic frequency of the harmonic trap) versus distance r 
from the centre of the trap (in units of a ), for two values of the reduced 
temperature T/T c ( T/T c = 0.1 and 0.5). The dashed curves and the dots 
are for an ideal Fermi gas and for an ideal Boltzmann gas (at T/T c = 0.1) 
in the absence of bosons, respectively. The full curves are the results of the 
numerical solution of the full set of eqns (1) - (7). The values of the parameters 
for the various panels are as follows: (a) / = 0.57p, Nj = 10 3 and Nb = 10 6 ; 
(b) / = 0.57.g and N f = N b = 10 4 ; (c) / = 2g, N f = 10 3 and N b = 10 6 ; (d) 
/ = 2g and Nf — Nb — 10 4 . The values of the BEC transition temperature 
T c and of the Fermi temperature T/ are also indicated. 

• Figure 2. Fermion density n/(0) at the centre of the trap (in units of Oq 3 ) 
versus reduced temperature ksT ' /Tilj. The results of the simplified model in 
eqn (8) (long-dashed curves) are compared with those for an ideal Fermi gas 
(short-dashed curves) . The arrows mark the BEC transition temperature and 
the Fermi temperature of the mixture. The inset in panel (a) shows the tem- 
perature derivative of n/(0) (in reduced units) versus reduced temperature. 
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